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Abstract 



In the previous paper, new sets of parameters in place of the Michel parameters 
have been proposed to analyze the data on the muon decay /i'^ — > e'^i/e^l- Both 
{V — A) and {V + A) charged currents with the finite neutrino mass have been used. 
In the present paper, this parameterization is extended to the more general form, 
and the method of data analysis (least squares) is discussed to determine the rate 
of contribution from the {V + A) current. There is a simple form in which a set of 
parameters is related primitively with the physical quantities. It is shown that the 
Michel parameters are one of the other sets which are obtained from this simple form 
by rearranging one term in it. We derive the condition to get the equivalent information 
on the unknown physical quantities, when the data are analyzed by using these simple 
and rearranged forms separately. There is some possibility to get different results from 
these analyses, because the equivalent condition is very delicate and the QED radiative 
corrections should be treated carefully. We propose a consistent formula for the data 
analysis. It is useful to compare the value of the least squares for the simple form with 
the one for the prediction of the standard model, because the large difference is not 
expected, especially for the Majorana neutrino case. Although we proposed a method 
to discriminate the type of neutrino in the previous paper, it is shown that it is not 
correct. 
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§1. Introduction 



The normal muon decay has been studied as a tool with high statistics to determine the 
structure of the weak interaction. The purpose of this paper is to investigate the effect of 
the {V + A) current added to the standard model and to find a means to treat the QED 
radiative corrections in the data analysis based on the method of least squares. Both the 
Dirac and Majorana neutrino cases are examined. 

Recently, the TWIST group^^ has reported their precise experimental data and analyzed 
them by using the helicity preserving four fermion weak interaction with {S ± P), {V ± A) 
and T forms. Their expression based on the Michel parameters for the e"*" energy spectrum 
is as follows; 



oc 



dx d cos 9 



Af{x) + P^coseV{x)j, (1-1) 



where 



M{x)^6x^\{1-x) + 1{Ax-3)pm\, (1-2) 
V{x) = 2 Cm [(1 -x) + ^{Ax- 3) 5m] ■ (1-3) 

Here x is defined as a; = E/W, where E is the energy of the emitted positron and W = 
{w?^ + mg ) / (2 m^) , and and mg are the muon and positron masses, respectively. The 
angle 9 is the direction of emitted e+ with respect to the muon polarization vector at 
the instant of /i+ decay. In the above expressions, we do not include terms proportional 
to m-g and neutrino masses, and also QED radiative corrections in order to simplify our 
explanation. 

The standard model predicts pM = = 0.75 and = 1 for these Michel parameters. 
The traditional way of experimental data analysis has been to determine the deviations from 
these predicted values. The new experimental center values reported by the TWIST group^^ 
arc pm = 0.75080, Sm = 0.74964 and 0.9960 < P^^m < 6/ < 1-0040. The QED radiative 
corrections are taken into account in their analysis. As you see, these deviations are small, 
so that it is preferable to determine them directly. 

In our previous paper^^ which is referred to as the paper I hereafter, we proposed to use 
new parameters that is suitable for investigating these deviations. We showed that various 
parameterizations are allowed by adopting the possible choices of the normalization factor 
for the isotropic part of energy spectrum. Among them, we mainly discussed the specific 
one which is related with the Michel parameters directly. This will be referred to as the 
Michel parameterization. We assumed that the weak interaction Hamiltonian consists of 
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both (y — A) and (V + A) charged currents and that the neutrino has the finite mass. 
Therefore, we have two kinds of lepton mixing matrices and three weak coupUng constants 
which represent the rate of mixture of {V + A) current. These unknown quantities will be 
referred to as "the weak coupling constants" in short. 

In this paper, the more general parameterization is investigated within the frame of 
the same Hamiltonian. In the simple form, a set of parameters is expressed primitively 
in terms of the combination of the weak coupling constants. Meanwhile, it is shown that 
the specific form related with the Michel parameterization is one of many forms which are 
derived from the above simple form by rearranging a term in it. Sets of parameters in 
these rearranged forms are related with the combinations of the weak coupling constants in 
somewhat complicated manners. Of course, these different sets of parameters should offer 
the same information on the weak coupling constants. In §3, we investigate the condition 
to get the same information mentioned above, when the experimental data are analyzed by 
using the method of the least squares for these different sets individually. It is pointed out 
that this condition is very delicate and we should be careful in treating the QED radiative 
correction in the data analysis. 

In the paper I, we also proposed a method of discriminating between the Dirac and 
Majorana neutrino experimentally by using the method of least squares for the e"*" energy 
spectrum. However, we shall show in §2 of the present paper that this proposal is incorrect. 
This discrimination is not easy in the muon decay. It will be discussed in §4. 

In §2, we present the comprehensive discussions on the general form of parameterization. 
The present experimental limits are listed for new parameters. In §3, we discuss the condition 
to get the information on the unknown weak coupling constants and the method of taking 
the QED radiative corrections into account appropriately. We propose the consistent formula 
for the method of the least squares in the data analysis. In §4, we present some comments. 
Appendix A contains expressions for the polarization of the emitted positron. 

§2. Parameterization of the decay spectrum 

We assume the following form of effective weak interaction Hamiltonian for the /x"*" decay, 

Hwix) = ^ [jlj^^ + XJI^Xr + ^IrA + >^31lA] + H.c. , (2-1) 

where Gp is the Fermi coupling constant. Weak coupling constants (A, 77 and k) represent 
the rate of mixture of the {V + A) current for the combination of the left(right)-handed 
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charged leptonic currents jiL{R) These currents are defined as 

2n 

JiLaix) = - 75)f/,,iV,(x), (2-2) 

i=i 

2n 

Jma{x) = J2Ei{xha{l + j5)VejN,{x), (2-3) 
i=i 

for the case of the n generations. Here Ugj and Vgj are the left- and right- handed lepton 
mixing matrices, and Ei and Nj represent, respectively, the mass eigenstates of charged 
leptons and neutrinos. Throughout this paper, neutrinos are assumed to have finite masses. 
The decay spectrum of in the rest frame of polarized /i"*" is defined as 

= rwA [Afix) + P, cos e Vix)] , (24) 
where the sum over the spin of has been taken and 

12 TT-^ 

The isotropic and anisotropic parts of energy spectrum obtained from the leptonic Hamil- 
tonian in Eq. (1241) are expressed as follows:^-*' 

A/'(x) = [a+ni{x) + {k+c + £mk+m)n2{x) + A c/^ ^3(2;)] , (2-6) 

Vix) = [a^pi{x) + {k^c + emk-m)P2{x)] , (2-7) 

where the decay formulae for the Dirac and Majorana neutrinos are obtained by setting 
Em = and Em = respectively. 

The constant A is introduced in Eqs. (1241) . (I2-6P and (12 -yp to simplify the coefficient 
of the prediction by the standard model in A/'(x). Its explicit form will be given later, see 
Eq. (12-231) . This A will be called a normalization factor according to the paper I. 



In order to see the physical meaning of weak coupHng constants, let us consider a typical example 
of the gauge theory, that is, the SU{2)l x SU{2)r x U{1) model with left- and right-handed weak gauge 
bosons, Wl and Wr. Weak coupling constants in Eq. (|2-ip are related to the physical quantities: 

X ^ {\c + t&v? C) , K — -q ^ {— taxi C,) . 

Here Ac = (Afi/M2)^, where Mi and M2 are masses of the mass-eigenstate gauge bosons which are expressed 
in terms of the weak eigenstate gauge bosons Wl and Wr with the mixing angle ^. For example, see Appendix 
A of the paper I. 
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The x-dependent parts are defined as follows: 



ni{x) 
n2{x) 

P2{x) 



(2-8) 
(2-9) 
(2-10) 
(2-11) 
(2-12) 



12a;pX (1 — x), 
6xpXo (1 — x), 



xli-l + 2x-rl), 
12xl{l-x), 



where 
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9.7-10 



-3 





4.7-10 



,-5 



(2-13) 



W 



The first terms ni{x) and pi{x) in Af{x) and V{x), respectively, are predictions from 
the standard model. They are named the standard functions in this paper. Meanwhile, the 
others, n2{x), n3(x) andp2(a;), are called the deviation functions. In these functions, all terms 
proportional to the neutrino mass are neglected, because of the smallness of {mu/m^) < 
9 - 10~^. Here, rriu stands for a typical mass scale of emitted neutrinos and is taken to be 



Coefficients a±, k±c, k±m and dr in jV{x) and V{x) are constants. They are obtained 
by summing up some products of the weak coupling constants (A, r] and k) and the lepton 
mixing matrices {Ugj and Vej) only over the emitted neutrinos. In the Dirac neutrino case, 
they are defined by 



Here it is assumed that all neutrinos can be emitted in the n decay, and we have used the 
relation 



from the unitarity condition, because j in the sum runs over all n neutrinos. 

*' The spectrum functions J\f{x) and V{x) in Eqs. (|2-6p and (|2-7p are precise except in the very tiny 
range of order (m^/m^)^ < O(10~^^) from the maximum of x, Xmax, which is expressed as Xmax = 1 ^ 
[{rrij + ruk)'^ / {2 m^W)] by using masses of two emitted neutrinos {nij and mfc). This is because two kind 
of kinematical factors come from the phase space integral over emitted neutrinos and they can be set unity 
in the almost entire range of x. In other words, they show significant x-dependence only in the extremely 
narrow range of order {m^/m^Y near Xmax- It should be noticed that M{x) and V{x) become zero sharply 
at ajmax due to these kinematical factors. For details, see §2.1 of the paper I. 






(2-14) 



(2-15) 
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By contrast, in the Majorana neutrino case, we assume that there exist additional n 
heavy neutrinos which are not emitted in the decay. Then, these coefficients are given as 
follows: 



a±=[{l-u;^){l-uj:^)±X^%^vj;^], (2-16) 

k±c= (0 [K^i-Tre')v;:^±r]^-e^i-u;;^)] , (2-i7) 

k±m= (l][^^\w^\^±v'\wi;^\'] , (2-18) 



(2-19) 



Here, m^^, f^^, Wefi and Wefih are all small quantities which stand for the extent of deviations 
from the unitarity condition due to the existence of heavy neutrinos. 

s;|f/,,f = i-M^^ E;|v^,,f = w', (2-20) 

V^j = vJ^, T^l Vek U^k = (2-21) 

where the primed sum is taken over only n light neutrinos out of 2n neutrinos. Their orders 
of magnitudes are ~ Vi'^ ~ 0{{m^D/m^jiY) and ~ ~ Oimyn/myR), if the 

seesaw mechanism is assumed|*l] Here, and m^R are, respectively, representatives of 
Dirac type and right-handed Majorana type masses in the neutrino mass matrix. 

2.1. Isotropic part of the spectrum: J\f{x) 

Let us consider the isotropic part N'{x). For the aim to survey the deviation from the 
standard model, it is suitable to examine M{x) by treating the standard function ni{x) as 
a base of analysis. In order to see the relation with the Michel parameter introduced in 
Eq. fll-2p . we rearrange ni{x) in N'{x) as follows: 

A/'(x) = ^|ni(x)[a+ + s {k+c + Bmk+m) +temXdr] 

+ [n2{x) -sni{x)] {k+c + Bmk+m) + [nsix) - tni{x)]emXdr\, (2-22) 



where s and t are some arbitrary numbers. Here the normalization factor A is set to the 
following Ast to simplify the coefficient of ni{x); 

Ast = a+ + s {k+c + £mk+m) + tem\dr > 0. (2-23) 
It is natural to restrict s and t to the values satisfying the condition Ast > ^■ 



For details, see Appendix A and §2.2 of the paper I as an example. 
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Thus, the isotropic part is denoted as Mstix) which takes the following form: 

Nst{,x) =ni{x) + [n2{x) - sni{x)]p^''^ + [n^^^x) - tni{x)]r]^''\ (2-24) 
where two parameters p*-**^ and r^*^'^*) are defined as 

^(st) ^ fc+c + gmfc+m ^ ^(st) ^ £mAc[r^ (2-25) 

It is worthwhile to note that positive within the frame of Hamiltonian in Eq. (12-ip . 

because its numerator is positive, as seen from Eqs. (12-141) and fl2-16p - (12- 181) . 

Two combinations [n2{x) — sni(a;)] and [n-i{x) — tni{xy\ in Eq. fl2-24p play roles of devi- 
ation functions for Mst{x)- In the paper I, there were some misleading discussions; that is, 
deviation functions presented there are not linearly independent as they are, and all of their 
coefficients cannot be settled independently *i] 

The simple choice of s and t is (s, t) = (0, 0). Then, the isotropic part is expressed as 
follows: 

Koix) = n^{x) + n^ix) p(°°) + n^ix) r^^^^). (2-26) 

This is nothing but the original expression given in Eq. fl2-6p with A = Aqq = a+. It should 
be noticed that A 7^ 1 in principle within the frame of our Hamiltonian. That is, if the 
right-handed charged weak current or the existence of heavy Majorana neutrino is assumed, 
we have the following expressions from Eqs. fl2-14p and (I2-I6p respectively: 

Aoo = (1 + A^) > 1 for the Dirac neutrino case, (2-27) 

^00 ~ (1 — — ^V^) ^ Majorana neutrino case. (2-28) 

The parameters p^^^^ and rj^'^^^ are related with the weak coupling constants. For the 
Dirac neutrino case, they are 

pm ^l{n^ + n^)>0, (2-29) 

r/(°°) = 0. (2-30) 

For the Majorana neutrino case, they are 

p(00)~i[«:2 ^_2^|_|2^) ^_2^|^|2^] (2-31) 

r7(°°) ^ ^ARe(wv* WV^). (2-32) 



*' In the paper I, a normalization factor is denoted by Ane- A shortcut to reconstruct hnearly inde- 
pendent deviation functions is to set n = £. Presentations there should be corrected as A„i — > A„n and 
Pm, 0. The old A„„ corresponds to A2n0 in the present paper. 



8 



In these expressions, only the lowest order terms are kept by taking ^oo ~ 1- Note that we 
can get no direct information on A from the isotropic spectrum ^/st{x) in the Dirac neutrino 
case. By contrast, in the Majorana neutrino case, the parameter proportional to A. 

But the order of magnitude of both p(°°) and themselves seems to be very small, as 
seen from Eqs. (I2^ and ([232]). 

Concerning the relation with the Michel parameter pM, the relevant term in Eq. fll-2p 
can be reproduced from the spectrum Mst{x) of Eq. fl2-24p by taking (s, t) = (2, 0). That 
is, we have the deviation function [n2(x) — 2ni(x)] and its associated parameter p^^^^ as 
follows; 



n2 (x) — 2 ni (x) ^ 2x^ (3 — 4 x) , 




(2-33) 



It should be noted that the behavior of [n2(x) — 2rii(x)] in Af2o{x) is different from n2(x) ~ 
12 x^ (1 — x) in Mqq^x). This will be discussed at the end of §4. 

The full expression of the Michel parameterization contains another combination 71,3 (x) tjm 
which is omitted in Eq. fll-2p because n^lx) is small due to a factor xq, as seen from Eqs. fl2-10p 
and fl2-13p .^^ There is a corresponding combination n3{x)r]^'^^^ in J\f2Q{x). These two pa- 
rameters, r]M and r]^'^^\ originate from different theoretical models for the weak interaction, 
as will be discussed in §4. But they have the same experimental values phenomenologically. 
We list here the experimental results reported by the Particle Data Group:^^ 

p(20) = _(6±7).io-4, (2-34) 
7^(2 0) = (^li 24) • 10^1 (2-35) 

These parameters are related with the weak coupling constants through the same ex- 
pressions as Eqs. fl2-29p - fl2-32p within the lowest order approximation. This is because the 
normalization factor A20 takes the following form: 

•A20 = (1 + A^ + + r/^) > 1 for the Dirac neutrino case, (2-36) 
~ [1 - - + (iV^ + + T]^ {v;^ + )] ^ 1 

for the Majorana neutrino case. (2-37) 

By the way, the spectrum Afst{x) with s 7^ and/or t 7^ is derived from Afoo{x) 
by rearranging the standard function ni{x) in it. Therefore, they should offer the same 
information on the weak coupling constants (A, r] and k). Indeed, this situation is expressed 
formally as the following identity: 

AtMt(x) = AoAroo(x), (2-38) 
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as seen from Eqs. f l2-23l) and (12 •241) . But the normalization factors Ast and ^oo themselves 
do not appear in the data analysis, and furthermore the deviation functions are different in 
the spectrum functions A/'sj(x) and f/oolx). We shall examine in §3 how to obtain the same 
information mentioned above, when the experimental data are analyzed by using different 
spectrum functions. For this purpose, the notation Afstix) will be used to indicate the 
spectrum with s 7^ and/or t ^ hereafter. 

There are some useful interrelations between parameters in jVoo{x) and those in jVstix). 
The following two are obtained directly from definitions in Eqs. (12-231) and (I2-25P : 

P^''^ = (i_,p(.*)_t^(.t))P^'*^ (2-39) 
V^"'^ = T. A TTTT^^'*^- (2-40) 

By using these relations, numerical values of p(°°) and r^(°°) can be estimated from the 
experimental results on p*^^°^ and rj'^'^^\ There is another identitjl*^ 

(l + sp(oo)+^^(oo)^) (i_,p(^t)_t^(-t)) =1, (241) 
from which we can derive the inverse relations to express p*^'^*-' and r/*^'^*-' in terms of p^'^'^) and 

^(0 0)^ 

It can be shown by using these identities that the relation in Eq. (I2-38P can be expressed 
as follows: 

-^-(-) = (i + ,,(oo)+,,(oo)) ^oo(x). (2.42) 

This implies that the spectrum function J\fst{x) as well as parameters {p^^^\ t^C'^*)) can be 
obtained from knowledge about Mqq^x) and (p'^°°\ 77^°°)), and vice versa. 

2.2. Anisotropic part of the spectrum: V{x) 

Next, let us consider the anisotropic part V{x) in Eq. {\2-7\i . We take the standard function 
Pi{x) as a base of rearrangement. In order to see the relation with Michel parameterization, 
we define a common factor by using the coefficient of ^2(2;); 

Bu = a_ + u{k_c + £mk-m), (2-43) 



We can derive the following relations from definitions in Eqs. (I2-23P and ()2-25p : 
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where u is some arbitrary number o Then, the anisotropic spectrum is denoted as Vstu{x)'- 
Vstuix) = e^^*") + [p2ix)-up,ix)] , (244) 

where parameters are defined as 

For the simple choice (s, t, u) = (0, 0, 0), we have ^00 = *^+ and Bq = a^. Then, the 
anisotropic spectrum is expressed as 

Pooo(x) =e^°°°) [pi{x) +P2(x)5(°)] . (2-46) 

This Pooo(a;) is identical with Eq. ([231) by taking = (a_/a+). 

The parameters and S^^^ are related with the weak coupling constants. For the 

Dirac neutrino case, they are 

e(°°°) = |i^, 5(°)-l(«:^-r^^). (2-47) 

For the Majorana neutrino case, they are 

^(000) ^ 1 _ 2X'v;%' ^ 1, ^ 1 [k' + - r^' {v;' + \w^,\')] . (2-48) 

Here, only the leading terms are given for the deviation from the standard model. 

On the other hand, the Michel parameterization in Eq. (II -SI) can be reproduced from 
Vstu{x) by choosing (s, t, u) = (2, 0, 6); that is, the following correspondences are obtained: 

P2{x) - 6p^{x) ^ 6x\3 - Ax), e^206) ^ ^^^^ ^(6) = ^ (^^ - 6m^ • (2-49) 

The experimental results reported by the Particle Data Group^-* are as follows: 

1^(206) ^ ^ QQ27 ± 0.0079 ± 0.0030, (2-50) 
5(6) = (1.1±2.7) ■ 10"^ (2-51) 

Here stands for the longitudinal polarization of the muon introduced in Eq. (Il-ip . The 
common parameter ^(^o^) is related with the weak coupling constants: 

^(2 6) = ^ ~ ^ '''f ~ for the Dirac neutrino case, (2-52) 
I + + K'^ + r]^ 

^ 1 + 2k\v;;' + w^') - A^\vf + tzv^^) ^ 1 

for the Majorana neutrino case. (2-53) 



*^ In the paper I, a common factor is denoted by Bni- Presentations there should be corrected as 
Bni, — > Bnn and 5m — > 0. Note that corresponds to B2n in the present paper. 
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The parameter 6^^^ is expressed in the same form as S^^^ in Eqs. (IMTj) or (I23HD if only the 
leading terms are kept for the deviation from the standard model. 
Parameters in Vstu{x) and ^000(2;) satisfy the following identities: 

f(ooo) ^ [1-^(5^"^] 

^ _t^(.t)]^ ' ^^^^^ 

The inverse relations are obtained by using the following identity. 

=1. (2-56) 

All these relations are derived from the definitions in Eqs. fl2-43p and fl2-45p . Also we can 
confirm the following relation, 

^^^'^^^^ = (l + ,p(oo) +^^(00)) ^000(3:). (2-57) 

Note that this relation is independent of u introduced in Eq. (12 -431) . 

Finally, it is useful to note that we have k = t] ii the SU{2)l x SU{2)r x U{1) model is 
assumed, and that the 5'-"'* parameter becomes simpler: 

= for the Dirac neutrino case, (2-58) 

^ ^ [{-,'--e') + (i^p - i^r)] « 1 

for the Majorana neutrino case. (2-59) 

§3. Method of least squares 

Let us find the condition to get the same results for the weak coupling constants (A, r] and 
k), when different spectrum functions are adopted in the data analysis. Also, the method is 
examined to include the QED radiative corrections. We shall first concentrate our discussion 
on the isotropic part Afst{x). It is easy to extend our treatment to the full spectrum where 
the anisotropic part Vstu{x) is taken into consideration. 

3.1. Analysis of the isotropic part of the spectrum: A/'(x) 

We shall use the method of least squares in the data analysis. The QED radiative 
correction is not included for a while. In the case where JVstix) is used, the unknown 
parameters (p'-'**-' and ry(**)) are settled as some values, at which the following x^st takes a 
minimum: 

1 r 1 2 

xlt = y^^— ^i^i) - CstJ^stixi) . (3-1) 
^-^ cr,- L J 

I) t 



12 



The summation over i runs over all measuring points Xj. The notation S{xi) stands for 
an experimental datum at Xi, and cij is its experimental error. The global normalization 
constant Cgt is introduced to adjust the theoretical values to the experimental data, so that 
the minimum point of x^st sought under the variation of Cgt as well as the parameters. 

By requiring that x^st takes a minimum, a set of analytical solutions is obtained for 
parameters (cst, Cst p^^^^ and CgtV^^^^)- They are known as the Cramers' formula for the 
system of linear equations, because these parameters appear linearly in CstJ^st{x)- If the 
spectrum function Afoo{x) is adopted and Xoo required to have a minimum, parameters 
(coO) CooP^°°^ and Coo^''°°0 are expressed by another set of analytical solutions. These two 
sets of solutions indicate that two global normalization constants Cst and Cqo satisfy the 
relation: 



It can be confirmed further that, under this relation, solutions for (p'-'^*-', p^^^^) and (r]^'^'^\ 7^(00)^ 
are consistent with their interrelations in Eqs. fl2-39p and fl2-40p . respectively. Therefore, the 
relation in Eq. (13 ■21) will be called the equivalent condition hereafter. 

If we combine Eq. (13-21) with Eq. (12-421) . the following equality is obtained: 



This means that the x^- values are the same for different spectrum functions A/'st(x) and 
Afoo{x); that is. 



In summary, the parameters (p*^°°^ and r]^^^^) are settled experimentally, when the Xoo" 
value becomes a minimum. The consistency between independent data analyses with use of 
Mooix) and Afst{x) is guaranteed by the equivalent condition for the global normalization 
constants (cqo and Cst)- Then, x^st becomes to be equal to Xoo) when their parameters 
satisfy interrelations in Eqs. (12-391) and (l2-40p . But, we should be aware that the equivalent 
condition is on the delicate balance of the global normalization constants, cqo and Cst, because 
their difference is very slight due to smallness of p(°°) and r]^^^\ 

Now let us examine how to take account of the QED radiative correction in the analysis 
of the spectrum. As the first step, we shall consider the case where the data analysis is 
performed by assuming the standard model. Then, we try to find a minimum value of the 
following X^^: 




(3-2) 



CstMstix) = cooAfoo{x). 



(3-3) 





(3-4) 




(3-5) 
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where f{x) stands for the QED radiative correction associated with the standard function 
ni{x) in Eq. fl2-8p .^ri Note that the unknown parameter in X^^ is only the global normal- 
ization constant Cgm- 

If the effect due to the weak coupling constants (A, r] and k) is considered, then the above 
XIj^ is modified. The standard function ni{x) is replaced by the spectrum function A/oo(a;) or 
J\fst{x) in Eq. fl2-24p . but the QED radiative correction f{x) is received no influence because 
of the consistency of the approximation. Thus, it is appropriate to introduce the following 
form instead of X^^: 

^It = E A [^(^0 - Cst^^st{x^) - C« f{Xi)f . (3-6) 

Here it is understood that the notation {s, t) includes the case of (s = 0, t = 0). Two new 
parameters Cst and Cr are different from Csm in -'^sm) because J\fst{x) is present in place of 
ni{x). Note that this is independent of s and t. This can be confirmed by comparing 
two sets of analytical solutions for parameters (cqo, p^^^\ 7^(00)^ c^) and (c^t, p^'^*\ ^^C**)^ c^). 
Each of these sets is obtained by requiring that Xqq or X^^ has its minimum. 

The equivalent condition in Eq. fl3-2p is derived again for these new global normalization 
constants, Cqo and Cgt, introduced in Eq. fl3-6p . Then, it can be proved under this equivalent 
condition that the x^-values are the same for different spectrum functions Afoo{x) and Afst{x) 
whose parameters satisfy interrelations in Eq. fl2-39l) for {p^^^\ p^^^^) and those in Eq. fl2-40p 
for (77(0°), r/(^*)): 

Xlt = ^OQ- (3-7) 

It is worthwhile to make the following three comments. First we shall examine whether 
the special choice of the spectrum function is preferable in the actual numerical analysis. 
For this purpose, it is useful to estimate the x-dependence of the function A{x) defined by 

A{x)=S{x)-c,m[niix)+fix)], (3-8) 

where Csm is fixed such that X"^^ in Eq. fl3-5l) takes a minimum. If this A{x) shows any x- 
dependence clearly, then we may choose such a value of s that A{x) is roughly proportional 
to the deviation function [n2{x) — sni{x)] in Eq. (12-241) . However, it is imaginable that A{x) 
does not show any clear x-dependence, because of experimental errors. If this is the case, 
then it may be preferable to adopt Afoo{x) in Eq. ( ]2-26p . because parameters are related with 
weak coupling constants in simpler forms. As a conclusion, we propose to use Xqq in the 
actual data analysis for its simplicity. 



r 1 

Relation between our f{x) and F{x) by Arbuzov^^ is f{x) — 1 + (mg/m^)^ — M2o{x). 
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The next comment is that the precise determination of the Xgj„-value itself defined in 
Eq. fl3-5p is interesting. This is because the large deviation from the standard model cannot 
be expected, especially for the Majorana neutrino case. This will be discussed at the last 
paragraph in §4. 

The final comment is that, in contrast to X^^ in Eq. fl3-6l) . the following definition 

^It = E ^k(^*) - [^-*(^*) + /(^^)] (3-9) 
^—^ cr. I 

is not appropriate theoretically, because it leads to the inequality Y'j^ ^ Fqq. 

3.2. Analysis of the full spectrum: V{x) 

In the extended form of the parameterization, the full spectrum in Eq. (l2-4p is expressed 
as follows: 

Vstuix, e) = [Afstix)+P^ COSeVstuix)^, (3-10) 

where JVstix) and Vstu{x) are given by Eqs. fl2-24p and fl2-44p . respectively. The method of 
least squares can be applied similarly to the isotropic part Afst{x). 

We summarize the essential points for the case with no radiative correction. The new 
X^stu defined as follows: 

1 r -1 2 

xltu = 2^— £{xi, Oj) - dstVstu{xi, 9j) , (3-11) 

where Xi and Qj are a set of observed quantities at one measuring point. We have the 
analytical solutions for two new parameters (^(**") and 5^"-') in addition to three old ones 
[dsti P^^^^ and ry(^*)) by requiring that Xstu takes a minimum. The corresponding solutions 
are also obtained by treating Xooo- should be noted that the global normalization constant 
d^t here depends only on s and t, because it is settled as a coefficient for Mstix) of Vstu{x, 6) 
in the method of least squares. Then, it can be verified not only that the global normalization 
constants satisfy the similar equivalent condition to the one in Eq. f l3-2l) . but also that other 
four parameters are consistent with interrelations in Eqs. (12-391) . fl2-40p . (12 -540 and (l2-55p . 
After all, due to this equivalent condition together with relations in Eqs. (12-420 and (l2-57p . 
we have the identity dst1^stu{x, 6) = dQQVQQQ^x, 0) and subsequently the equality 

xltu = Xooo- (3-12) 

In the case where the radiative QED effect is taken into consideration, we modify x\tu 

and define the following Z'].^^ which satisfies both the consistency conditions for parameters 
2 _ 72 . 

stu ~ ^000- 

Zlt^ = '^— £{xi, 9j) - dsiVstuixi, 9j) - dj,F{xi, 9j) , (3-13) 



and the equality Z^^^ = Z"^ 



■ ■ 
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where a new parameter corresponds to in Eq. f l3-6p and 

F{x,e)= [fix) + P, cose g{x)]. (3-14) 

Here g{x) stands for the QED radiative correction associated with the anisotropic standard 
function pi{x) in Eq. mihM 



Corresponding to for the isotropic part of the standard model, the following Z 



2 



sm 



with a global normalization constant dg^i is defined for the full spectrum: 

^sm = 5^ ^ - ^^sm [DsmiXi, Oj) + 9,)] \ (3-15) 



I, J u 



where 

Am(x, 6) = [ni{x)+P^ cosOpiix)]. (3-16) 

It is interesting to compare minima of Z^^ and in order to know directly the extent of 
departure from the standard model. 

§4. Discussion 

Let us consider the possible method to determine whether the neutrino is of the Dirac 
or Majorana type. It is offered by observing the t]^^^^ parameter which is zero or nonzero 
depending on the Dirac or Majorana neutrino within the frame of the gauge theory, as 
seen from Eqs. (12-301) and (12-321) 1**^ Meanwhile, the Michel parameter t]m has been popular 
as a measure to show the deviation from the standard model. It corresponds to r]^'^^\ as 
mentioned in Eq. fl2-35l) . But, the r/^**^ term is defined for the Majorana neutrino case 
within the frame of gauge theory, while the tjm term comes from the interference between 
the {V ± A) and (S* ± P) (or T) forms even for the massless Dirac neutrino case.^-* Anyway, 
the observation of this r] parameter indicates the deviation from the standard model. 

It is well known that the detection of this rj parameter is very difficult experimentally. 
One of its reasons is that the contribution of the relevant deviation function n^lx) is small, 
because it is proportional to the small value of xq, as shown in Eq. (I2-I3p . We may avoid 
this weak point by considering the r-decay: 



r 



/i+ + z/^ + z77. (4-1) 



r 1 ~^ 

*' Relation between our g{x) and G{x) by Arbuzov*-' is g{x) = — 1 + (me/m^)^ G{x) — V2q{x). 
**) This difTerence is independent of the choice of the normahzation factor A. It should be noted on this 
point that there are some misleading discussions in the paper I. It is not correct to say that there is some 
difference between the Dirac and Majorana cases by choosing A. 
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All formula in the previous sections can be applied to the r-decay by the replacement of 
both (m^ rrir) and {rrie Tn^). Then, the value of Xq = me/W ~ 0.01 is shifted up to 
xo{t) = rrin/Wr — 0.12 where Wr = {rn^ + m'j^/2m-r- Meanwhile, the second reason for 
the difficulty to detect the r^^^*^ parameter is its smallness within the gauge theory, as seen 
from Eqs. fl2-30p and fl2-32p . Its rough estimate was discussed in §4.2 of the paper I. As a 
conclusion, the muon decay cannot be used to discriminate between the Majorana and Dirac 
neutrino cases in reality. 

Next, let us examine the order of magnitude of the normalization factor A. In the 
traditional model which has been used to analyze the experimental results, it has been 
assumed that A = = 1-^'' However, in our model for the Dirac neutrino, as we know 
^00 = (1 + '^^) > 1 from Eq. fl2-27p . we need some information on A^. We may say that 

< O(10~3) from the reported data in Eqs. fl2^ - fl2^ and fl2^ . In the Majorana 
neutrino case, there is no definite information at present, although it is imagined from 
Eq. (12-281) that the deviation from ^00 = 1 is very small. 

Finally, we would like to comment on the data of the Michel parameter p^. The mean 
value of p^^^-* obtained from pm is p^^^-* = —6 ■ 10~^, as shown in Eq. (12 ■341) . This mean 
value is negative, although it can be positive within experimental uncertainty. From the 
theoretical side, it is predicted to be positive within the frame of Hamiltonian in Eq. (12-ip . 
as seen from Eqs. (12-290 or (12-3ip and ( l2-39p . There is a possibility that this difference 
comes from some ambiguity in the data analysis. This is because the consistency of the 
data analysis depends delicately on the equivalent condition in Eq. (l3-2p and the means 
to treat the QED radiative correction, as mentioned in §3. Under these circumstances, it 
is of interest to compare experimental results for p*^^°^ and p'^°°\ which are determined by 
using Af2o{x) and A/oo(a^), respectively. These parameters should satisfy the interrelation 
in Eq. (12-390 and have the same signature. In this connection, the evaluation of Xqq and 
X20 is of importance, because they are equal theoretically. It is also interesting to compare 
them with X^^ in Eq. (13-50 . because the large deviation from the standard model cannot be 
expected, especially for the Majorana neutrino case. 

Appendix A 

Polarization of positron 

Since parameters are defined in somewhat different forms from the paper I, we shall list 
expressions for the polarization of an emitted positron. 
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The differential decay rate is expressed as follows; 



i + p,ix,e)-c , (A-i) 

ax a cos 6 2 L J 

where the vector Pe{x, 6) is a polarization vector of e^, and ( is the directional vector of 
the measurement of the e"*" spin polarization. The decay plane is defined by the momentum 
direction (p^) of e+ and the muon polarization vector {Pfj.). 

Three components of the e"*" spin polarization vector are defined as^^ 

P^ix, e) = P,ix, e)p\ + PMX, X -^m) X ^ p^^^^^ ^) Aii^L. (A.2) 

\{Pe X Pf,) X Pel iPeXP^l 

The explicit expressions of these Pl{x, 9), Pti{x, 9), and Pt2{x, 9) are presented in terms 
of parameters defined in §2 of the present paper. For simplicity, they are listed only for the 
simple form with both the normalization factor ^oo and the common factor Eq. Also, the 
radiative corrections are not included here.®^ 

A.l. Longitudinal polarization: Pl{x, 9) 

It is convenient to separate the isotropic and anisotropic distributions of e"*" with respect 
to the muon polarization vector P^, namely, 

D{x, 9) 

where the denominator D{x, 9) is defined from Eq. ( 13-lOp as follows: 

D{x, 9) = —Voooix, 9) = —[Afoo{x) + P^ cos9Voooix)]. (A-4) 

The isotropic part Q{x) and anisotropic part S{x) of the longitudinal polarization are, 
respectively, expressed as follows: 

g(^)=^{ooo) [g,(x)+g2(a;)5(°)], (A-5) 

S{x) = [si{x) + S2{x) p(°°) + S3{x) r/(°o)] . (A-6) 

where 

qi{x) = Xp{3 -2x - rl), (A-7) 

q2ix) = 12xp{l-x), (A-8) 

= (-X + - Xq), (A-9) 

S2{x) = 12x{l-x), (A-10) 

S3(x) = -2xo(l-x). (A-11) 

The parameters (^(°°°) and 5*^°^) in Q{x) are defined in Eq. (12-451) for the case of V{x), while 
the parameters and r/*^°°^) in S{x) are defined in Eq. (I2-25I) for the case oiJ\f{x). 
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A. 2. Transverse polarization within the decay plane: Pti{x, 9) 
The x-dependent part R{x) of Pti{x, 9) is defined as 



with 



R{x) = [n(x) (1 - 12p(°»)) + r2(i) !(<»"'] , (A-13) 



where 



ri(a;) = -,To (1 - x), (A-14) 
r^i^x) = -2{x-xl). (A-15) 

Note that the small quantity xq appears in ri (x) , which stands for the prediction from the 
standard model. Meanwhile, rj'^^^^ which indicates the existence of the Majorana neutrino is 
associated with the larger deviation function r2{x). 

A. 3. Transverse polarization perpendicular to the decay plane: Pt2{x, 9) 
The x-dependent part T{x) of Pt2{x, 9) is defined as 

P^ sm9T{x) 

^' - D{x,9) ' ^^'^^^ 

where 

T{x)=2x,{l-rl)7j'iT. (A.17) 
Here the new parameter rj^^^ is defined as follows 

viT = £m (^4^) Ini(l^* IZ^). (A-18) 

A non-zero value of T{x) imphes the existence of a non-zero Majorana CP violating 
phase in our model. There is no corresponding term in either the standard model or our 
model for the Dirac neutrino. 
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Errata added in September, 2009 

1. The inequahty > in Eq. (2-25) should be replaced by p^"*) > 0. 

2. Equation (2-36) in §2 should be changed from 77(^0) = (1 ± 24) • 10"^ to 

,y(20) = -0.12 ±0.21 or - (2.1±7.0± 1.0) • 10~^ for the Major ana neutrino case. (E-1) 

Here, the former rj^"^^^ — —0.12 ± 0.21 is reported by Derenzo^^ from his analysis of the 
energy spectrum. The theoretical expression he used is identical to our N2o{x) in Eq. (2-24). 
On the other hand, the latter is obtained by Danneberg et al.^^ through their restricted 
analysis of the transverse polarization Pti{x, 9 = n/2) of e"*". Their approximated theoretical 
expression is the same as ours under the assumption p^'^^\x) = 0, because they use the one- 
parameter fitting by assuming pm = 0.75 (i.e., p*^^°^ = 0). Our expression of Pti in the 
(s, t, u) = (2, 0, 6) case is expressed without any approximation as 

PM.,e = ./2)= ^,^^ ,1;^,) . (E.2) 
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where 

^ = 7r/2) = (3x - 2x^ - x^) + 2 (3x - 4^' + x^) p^'o^ + 6x0(1 - x)/^^'"^ (E4) 

Note that our expression of Pti in the (s, li) = (0, 0, 0) case is given in Eq. (A- 12) in 
Appendix A. 

The data 77 = (1 ± 24) • 10~^ that we cited in our paper is that reported by the Particle 
Data Group^^ as an average over the values of 77 obtained by different experiments. This 
citation is not appropriate, because some of them are derived from the data on Pti{x, 9) by 
using theoretical expressions different from our R^'^^\x) in Eq. (E-3). 
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